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Linear stability analysis of the Hall magnetorotational instability in a 
spherical domain 
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We investigate the stability of the Hall-MHD system and determine its importance for neutron stars at their birth, when they 
still consist of differentially rotating plasma permeated by extremely strong magnetic fields. We solve the linearised Hall- 
MHD equations in a spherical shell threaded by a homogeneous magnetic field. With the fluid/flow coupling and the Hall 
effect included, the magnetorotational instability and the Hall effect are both acting together Results differ for magnetic 
fields aligned with the rotation axis and anti-parallel magnetic fields. For a positive alignment of the magnetic field the 
instability grows on a rotational time-scale for any sufficiently large magnetic Reynolds number. Even the magnetic fields 
which are stable against the MRI due to the magnetic diffusion are now susceptible to the shear-Hall instability. In contrast, 
the negative alignment places strong restrictions on the growth and the magnitude of the fields, hindering the effectiveness 
of the Hall-MRI. While non-axisymmetric modes of the MRI can be suppressed by strong enough rotation, there is no 
such restriction when the Hall effect is present. The implications for the magnitude and the topology of the magnetic field 
of a young neutron star may be significant. 
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1 Introduction 

Neutron stars are usually formed as a result of the gravita- 
tional collapse of a star at least eight times more massive 
than the Sun. Degenerate electron pressure, diminishing as 
the most energetic electrons begin interacting with the pro- 
tons in the nuclei producing neutrons and neutrinos, is un- 
able to stop the gr avitation al collapse until the nuclear den- 
sities are reached (IShapiro & Teukolskv 1983 ). 

The transition from the collapsing stellar core to a sta- 
ble neutron star occurs through the protoneutron star stage . 
At the end of the collapse, the remnant stabilises into a hot, 
superdense, degenerate nuclear fluid confined beneath the 
bounce shock zone of the ensuing supernova (Page et al. 
I2OO6.) . The proto-neutron star (PNS) has initially a diameter 
of around 100 km which reduces to the canonical 20 km in 



one second dBurrows & Lattimer 1986). The final transfor- 



mation of the core remnant into a young neutron star lasts 
for additional one hundred seconds, during which it settl es 
in the state of cold catalysed matter JBonanno et al.ll2006h . 
Simulations of the core collapse show that almost 
any initial rotation profile, includi ng rigid rotation, re- 
sults in a differentially rotating PNS ( A kiyama et al.ll2003 ; 



ObergauUnger et al.l I2006t lOttet al.. .20061) . Depending on 
the initial spin of the progenitor core, the collapse may 
result in a strong differential rotation with a negative 
radial gradient. Such conditions are favourable for the 
occurrence of the mag netorotational instability (MRI) 



(lBalbus&Hawlevlll99lf) 



Because of the very short time-frame during which hy- 
drodynamical instabilities can take place in the PNS, the ex- 
ponential growth of the magnetic field provided by the MRI 
may be an essential ingredient in bridging the gap between 
the relatively weak progenitor magnetic fields and those ob- 
served in strongly magnetised neutron stars. It may also play 
an important rol e in the mechanism which powers super- 
nova exp losions ( Akivama et al. 2003; Ardelian et al. 2 0051: 
Burrows et al. 2007: ObergauUnger et al. 2006. 2009) . 



The magnitudes of the magnetic field during the PNS 
stage of up to lO'^^ G raise the question of how important the 
Hall effect may be for its evolution. So far, the theory of the 
interaction between the MRI and the Hall effect has mostly 
been flourishing in the context of weakly ionised accretion 
disks. 

The local analysis of the Hall-MHD in differentially 
rotati n g disks ( Balbus & TerquemlbOOlt lUrpin & Riidiger 



2005 



I disks ( Balbus & Terquemll2001t lUrpin & Riidiger 
Wardld 119991) revealed that the influence can be 
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both stabilising or destabilising, depending on a number 
of factors, including th e orien tation of the magnetic field. 
Riidiger & Kitchatinovl ( 120051) found similar results con- 
cerning the global instability in protostellar disks and deter- 
mined threshold magnetic fields for which the Hall induced 
instability would occur 

An interesting feature observed in the these works was 
that an instability occurs even when the conditions for the 
usual MRI are not fulfilled. For fields parallel to rotation 
axis, weaker than expected magnetic fields are suscepti- 
ble to the instability. In addition, even a positive differ- 
ential rotation gradient - stable against the MRI - can be 
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destabilised. This is an indication that there are new chan- 
nels through which the magnetic field taps into the energy 
of the shear when the Hall effect is present. This is even 
more obvious in models w ith pure shear without rotation 
dBeiarano et aLll201ll;lKunz.2008.). 



Kondic et alJ ( 120 llh examined the instability driven by 



the Hall effect itself, the shear-Hall instability (SHI), in the 
context of protoneutron stars. Like the MRI, it feeds off the 
shear energy and grows on the time-scale of the rotation 
period. They show that the SHI may, if the conditions are 
fulfilled, occur in the PNS, but they focus solely on solving 
the induction equation, thus leaving out the back-reaction of 
the field on the flow. 



Unlike ' Kondic et alj J2011h . this works includes the 
back-reaction leading to the interplay between the Hall ef- 
fect and MRI. Stability maps and growth rates of the global 
Hall-MRI instability follow from numerical calculations of 
the Hall-MHD equations in a spherical shell. These results 
are compared to the properties of a protoneutron star 

2 Model 

The coupling between the magnetic field and the flow re- 
sponsible for the MRI and the SHI is contained in the 
Navier-Stokes and the induction equations of the MHD sys- 
tem. In the incompressible and isothermal case, these read 



du 
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dB 
~dt 
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The variables u and B are, respectively, the velocity and the 
magnetic field, while p is the pressure. The Hall parameter 
Rb is the ratio of the off-diagonal to the transverse com- 
ponent of the conductivity tensor and scales linearly with 
the magnetic field. The viscosity is represented by v, mass 
density by p and the isotropic diagonal element of the elec- 
tric conductivity tensor by a. The unit vector e^ gives the 
orientation of the magnetic field. 

Equations ([T]i and (|2]i do not differ from those used in 
many other works dealing with weakly ionised plasmas. 
However, the physics behind the interactions of different 
plasma species is not the same. The two charged compo- 
nents considered implicitly in this model are fully ionised 
atoms and electrons. The Eq. (|2|i is what remains after all 
the terms of order m^jm^ are neglected, where m^i are the 
masses of electrons and ions. This can be interpreted as if 
the motion of the fluid is entirely due to the ions, while the 
current carriers are the much lighter electrons. If the mag- 
netic field is weak, or the collisional frequency between the 
species is large, Rb tends to zero and ordinary MHD is re- 
covered. The Hall effect is then a correction to the ordinary 



MHD which arises from a response of electrons to the mag- 
netic field. 

Because the aim of this work is the investigation of the 
conditions under which the instability occurs, the system 
can be linearised around a stable (or a very slowly evolving) 
state. We chose a spherical shell with radial boundaries at 
Ti = 0.7 and r^ — 1. The shear necessary for the operation 
of the MRI and the SHI is provided by a hydrodynamically 
stable differential rotation profile 

N21-1/2 



^{s) = f^o 1 + (s/so)' 



(3) 



where ilp is the angular velocity on the rotation axis, s 
the cylinder distance from it and sq = 0.5. The spherical 
components of the background velocity, Ur, uq, u^ are thus 
Uo{s) = (0, 0, sJl). A homogeneous magnetic field BtySz 
is imposed in the direction of the rotation axis. This choice 
prevents a coupling between the flow and the field in the 
absence of a perturbation. 

The domain of linear solutions is restricted to a spheri- 
cal shell of constant mass density interfacing to the current- 
free surroundings on both of its sides. This choice was made 
for the sake of numerical tractability of the problem, even 
though this is a considerable simplification for the inner 
boundary. 

It is now possible to write the linear form of the equa- 
tions. With the separation of the magnetic field into pertur- 
bations and background field, B = b + Bq, we can write 



^ = S'iV X b) X e, 
ot 



v/- 



(4) 



(5) 



itoxVxM + itxVxMo + PmAu, 

— = -RbV X [(V X b) X e,] - 
ot 

V X (V X b) + V X (mo X 5 + M X JBq), 

where all gradient terms are collected in V/, which will 
vanish as we are taking the curl of Eq. (|4) for the solution. 
The time is normalised to the diffusion time-scale t-q = 
R? JT], r\ = (? jiA.'Ka) and all the units of length to the radius 
of the outer boundary R of the star in Eqs. (JUi and (|5]l. The 
Lundquist number is the ratio between the diffusion time 
and the Alfven time, ta = ^/^TTpR/Ba), whence 

RBq 



S 



(6) 



The magnetic Prandtl number Pm is the ratio of viscosity 
to magnetic diffusivity and is set to unity in the computa- 
tions presented in this paper The magnetic Reynolds num- 
ber is Rm — R^rto/rj. Finally, Rb denotes the Hall param- 
eter fixed by the magnetic field Bq . As the density is taken 
to be constant, these parameters do not vary in the domain. 
The information about the conditions for the onset of 
the instability and the time-scale of its development is con- 
tained in the growth rate of the eigenvalue solution of the 
system (|4|-(|5l). A linearised version of the pseudo-spectral, 
spherical MHD code by Hollerbach (2000) was used to de- 
termine the growth rates and the stability criteria. 
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The results of a local instability analysis are useful 
to determine the numerical limitations of the global ap- 
proach. According to th e ideal MRI dispersion relation 
(IBalbus & Hawleyl Il998h . the length-scale of the fastest 
growing MRI mode, for a given magnetic field tends to zero 
as the shear becomes weaker An ideal numerical simula- 
tion would face a lower limit of the shear of showing insta- 
bility because of the non-vanishing numerical diffusivity. In 
our case, the presence of a physical diffusivity sets a lower 
limit on the shear (as well as on the magnetic field for given 
shear) to be MRI unstable. We can compare the growth rate 
of the instabilities with the decay rate of the diffusion for a 
given wavenumber k, Wd — —f] fc^, and thus derive a max- 
imum wavenumber fcmax which the unstable mode of the 
MRI should have to be excited at a given finite 77. From 
1 dfl 




:V 



r!n 



(7) 



(8) 



2V V2?/ 

for the assumed rotation profile and s — sq ~ 1/2. In our 
dimensionless units, this corresponds to fcmax = 13 at 
Rm = 1000 and fcmax = 23 at Rm ~ 3000. The resolution 
of 20 Chebyshev polynomials and 70 Legendre polynomials 
was enough to resolve all the corresponding length scales in 
the range of investigated values of Rm and S. 

3 Results 

The eigensolutions of the system (|4|i-(|5]) are divided in sep- 
arate classes defined by the transformation properties of 
the equations. The system is invariant under both rotations 
around e^ and reflection in the equatorial plane. The for- 
mer implies on e'™"^ behaviour in the azimuthal direction, 
while the latter distinguishes between antisymmetric field 
(accompanied by symmetric flow) and symmetric field (an- 
tisymmetric flow) behaviour in the (r, 6')-plane. The sym- 
metric field has the property that only its 6'-component re- 
verses sign under symmetry transformation, while the ra- 
dial and azimuthal components remain the same. The op- 
posite holds for the antisymmetric field. The antisymmet- 
ric and symmetric m = solutions are labeled with AO 
and SO, respectively, in this manuscript. The label AO thus 
means that 6r(r, 6*, 0) = —br{i',Tr — d,a(t)), bg{r,9,(t)) — 
be{r,'K — 9,a(f)), b^{r,9,(l)) = —b^{r,TT — 9,a(f)), while 
Ur{r, 9, (j)) — Ur{r, TT — 9, a(j)), ug{r, 9, (j)) = —ug{r, n — 
9, a4>), u^{r, 9, 0) = Uci,{r, n — 9, acj)) for arbitrary a. 

We note that there is a difference between symmetries 
of the linearised and the full induction equation. While the 
solutions of the linearised equation can be of either A or 
S symmetry ty pes, only the A t ype is possible in the non- 
linear case (see lHoUerbach & R iidiger 2002 for detailed ex- 
planation of the symmetries of the Hall induction equation). 
The absence of total equatorial symmetry has exciting con- 
sequences for the resulting topology of the magnetic field, 
as we shall explain in the final section. 



E 




Fig. 1 Pure MRI stability curves for axisymmetric pertur- 
bations {m = 0) which are symmetric (SO) and antisym- 
metric (AO) with respect to the equator 



AO; S=0 

SO; S=0 

AO; S=40 

SO; S=40 




Fig. 2 Hall-MRI stability curves for the SO and AO modes. 
The instability region is inside the curves. The influence of 
the Hall effect, for a fixed background field {S — const) 
is represented by /3 = Rb/\S\. Pure SHI corresponds to 
5* = case. 



AO; S=0 

SO; S=0 

AO; S=70 

SO; S=70 




Fig. 3 Hall-MRI stability curves for SO and AO modes. 
The instability region is inside the curves. The influence of 
the Hall effect, for a fixed background field {S — const) 
is represented by /3 = i?B/|'S'|. The solid and long-dashed 
curves are solely due to the SHI. 
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3.1 Stability maps 

The curves in Figs. [T] |2] and [3] show the marginal stabil- 
ity of the Hall-MRI system. In Figs.[2]and[3l /3 is the ratio 
between the Hall parameter and the absolute value of the 
Lundquist number, /3 — Rb/\S\. For non-quantising fields, 
it depends only on the sign of the magnetic field and the 
structural properties of the plasma. It can be thought of as 
a measure of importance of the Hall effect for a fixed back- 
ground field. 

The curves in Fig. [T] represent the pure MRI. Note the 
transition to stability for S* ^ 25 due to diffusion. Here, the 
influence of the field on the fluid is exerted through the 
Lorentz force, but there is no transformation of toroidal into 
poloidal magnetic field since there is no Hall effect. 

The limit of small S is the opposite case. For a very 
long Alfven time, there is only amplification of magnetic 
perturbations by the differential rotation. This case is repre- 
sented by the S — Q line in Figs.|2]and|3] In the framework 
of ordinary MHD, this would be a stable situation. As it 
can be seen, however, the Hall effect gives rise to the SHI. 
The shear feeds the energy into the toroidal field, and the 
Hall effect channels it back to p oloidal creating a po sitive 
feedback. This is the same as in iKondic et alJ (1201 ih . The 
consequence is that, due to the shear-Hall mechanism, there 
is an instability even in the range 5 = to 5 = 25. 

Outside of the diffusion zone, the stability is under 
larger influence of the interaction between Navier-Stokes 
and induction equation (positive /? in Figs. |2] and (3]) than the 
shear-Hall mechanism, with the critical Rm being closer to 
the pure MRI values. For a certain window of /3, the pres- 
ence of the momentum equation has a slightly stabilizing 
effect compared to pure SHI in the induction equation only. 

The results for B^ anti-parallel to the rotation axis are 
more dramatic. There is no pure SHI in this case, but the 
Hall effect still produces some important differences. The 
instability zone now also possesses an upper limit of Rm be- 
yond which the system becomes stable. On approaching the 
upper boundary, the length-scale of unstable perturbations 
increases, until it finally exceeds the domain size at the crit- 
ical lin e. A similar effect was noted in a local approximation 
by Balb us & Terqueml(l200lk . 

The fact that the opposite orientation of Bq and il is 
less stable than when the two vectors are aligned, is already 
known from accretion disk theory (Riidiger & Kitchatinov 
20051) . In this case, the Hall effect enhances the MRI mech- 
anism. However, for large enough |/3| at a given S or for 
small enough 5* at a given (3, the Hall effect suppresses the 
instability, rendering the system stable. Apparently, there is 
no way for the MRI to set in for substantially lower critical 
Rm with the help of the Hall effect. 

The comparison of Figs. |2] and |3] also shows how the 
"depth" (minimal unstable /3) of the instability zone in- 
creases, if the field grows to larger strengths (larger |5|). 

To summarise, if the field and the angular velocity vec- 
tor are anti-aligned, the Hall effect severely restricts the 
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Fig. 4 The relationship between the growth rate (nor- 
malised to r2o/(27r)) and the magnetic field. In this graph 
the ratio between the Hall parameter and the Lundquist 
number is kept constant. Therefore, changing S amounts 
to changing the strength of the Bq field. The magnetic 
Reynolds number is kept fixed. 



maximal magnetic field that can be MRI unstable. There 
is no such restriction for the positive orientation. This is a 
crucial result, calling for reassessment of any conclusions 
about the field patterns emerging in differentially rotating, 
strongly magnetised systems, derived assuming the validity 
of ordinary MHD approximation (for example, a body of 
work dealing with the MRI in PNS). 

3.2 Growth rates 

We varied the strength of the magnetic field but kept the 
ratio of Hall parameter to Lundquist number, j3, fixed. The 
resulting growth rates are shown in Fig. 2] 



Similar to the findings by iRiidiger & Kitchatinov 



(120051) . the Hall effect shifts the extrema of the graph with, 
again, evident stabilising influence for the positive orienta- 
tion (for S > 70, Rb > 0). Also the largest growth rate ob- 
tained for Hall-MRI is limited and smaller than the growth 
rates of the MRI for large enough S. For S < 100, the 
growth rates of the Hall-MRI is smaller for an anti-parallel 
orientation of the magnetic field and the rotation than for 
parallel orientation. Only for very strong anti-parallel mag- 
netic fields, the growth rates exceed the parallel case, but 
never exceeds the growth rate of the pure MRI for S ^ 1. 
The maximal growth rate was of the order of the inverse 
rotation period, the smalle st time-scale of the system, sim- 
ilar to what was found by iKondic et alJ (1201 ih for the SHI 
above. 

3.3 Field topology 

Figure |5] illustrates the topology of the toroidal compo- 
nent of the most unstable AO mode for Rm = 3000 in the 
situations when the pure SHI, or the pure MRI or both 
mechanisms responsible for the field amplification are at 
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(a) (b) (c) 

Fig. 5 Vertical cross-sections of the magnetic field (black-blue) and the velocity field (red-yellow), (a) Toroidal fields 
for Rb ^ lands ^ 0, pure SHI. (b) Toroidal fields for Rb ^ and S ^ 70, pure MRI. (c) Toroidal fields for Rb = 1 
and S = 70, Hall-MRI. 
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Fig. 6 Critical magnetic Reynolds number for the pure 
MRI and the Hall-MRI for the to = 1 mode for which 
the magnetic field is symmetric with respect to the equator 
The Lundquist number S is varied, while the Hall parameter 
Rb = 1 is fixed. 



work. The rotation axis and the background magnetic field 
are parallel here. The formation of zones which grow in 
number with the increase of the magnetic R eynolds num- 



ber (he re Rm = 3000) is akeady known from lKondic et al. 



(120111) . Another interesting feature, evident from comparing 
Figs.|5t and|5]3, is that the corresponding sections of the two 
wound-up toroidal field plots have the opposite signs. This 
effect is a consequence of the entirely different manner in 
which the MRI and the SHI generate the radial field compo- 
nent (which is then wound up into the toroidal field by the 
shear). Both mechanisms acting together will interfere and 
decrease the efficiency of the amplification. 

Additionally, since the Lorentz force is invariant under 
the transformation Bq — > —Bq, while the Hall term changes 
sign, after the transformation the Hall term will produce the 
toroidal field in the opposite direction. The stability limit is 
increased correspondingly for positive /3 in Fig. [3] as com- 
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Fig. 7 Critical magnetic Reynolds number as in Fig. |6]but 
for the m — 1 mode which is anti-symmetric with respect 
to the equator. 



pared to the (3^0 stability limit. This is why both processes 
interfere constructively in case of Bq < (anti-parallel to 
the rotation axis) and small enough /3. The stability limit 
goes below the /3 = limit in the left-hand side of Fig. [3] 
These res ults are comparable to the local Hnear stability 
analysis of lBalbus & TerquemI ( 1200 lb . 



4 Non-axisymmetry 

The marginal stability of an to = 1 perturbation was also 
studied in this setup. The critical magnetic Reynolds num- 
bers for pure MRI and the Hall-MRI are shown in Figs. |6] 
and |7] for the SI and Al modes, respectively, as a function 
of S. A typical feature for the stability of non-axisymmetric 
perturbations in differentially rotating system is the pres- 
ence of an upper limit for th e Reynolds number (see also 
Kitchatinov & Rudigen2010l for similar behaviour in disks) . 
At 5' = 40, for example, there is an instability window for 
SI between 1130 < Rm < 5410. There is no instability 
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Fig. 8 Surface map of the radial component of the mag- 
netic field for an the case of an unstable Al mode at 
Rm = 5000, S = 35, and /3 = 1. Light areas represent hr > 
0, while dark areas are br < 0. 



for S < 33, while the system is still unstable against ni — 
perturbations. For small fields, the MRI only produces ax- 
isymmetric field geometries. The Al mode is unstable for 
1030 < Rm < 5130 at S* = 40 and entirely stable against 
pure MRI for S < 30.5. 

The situation is again changed by the presence of the 
Hall effect. Figure |6] shows the critical Rm for Rb = 1. 
There is practically no longer any dependence on S, with 
critical Rm always being near 1000. The critical values are 
usually slightly above the marginal Rm for the m = 
modes, but can be lower for large Lundquist numbers S (e.g. 
compare 5 = 70 in Fig.|6]with f3 — 0.014 in Fig.|3]l. There 
is obviously more non-axisymmetry possible in this differ- 
entially rotating system if the Hall effect is present. 

A surface plot of the magnetic field pattern emerging 
from an unstable Al mode is shown in Fig. |8] Like for the 
axisymmetric modes, the eigenmode is concentrated near 
the equatorial plane. There is no indication for distinct mag- 
netic poles as one would expect for a tilted dipole. The 
TO = 1 instability will probably not directly serve as an ex- 
planation for the non-axisymmetric magnetic fields of pul- 
sars. 



5 Discussion 

The physical model of PNS can be related to the stability 
map in Fig. [T]using the code developed by Potekhin (2008) 
for calculating the electrical conductivity tensor of dense 
plasmas, by taking into account the electron-ion interaction. 
However, at T > 10^" K, the plasma of the protoneutron 
star envelope contains a distribution o f fully ionised atoms , 
free neutrons, protons and electrons (IHaensel et al.ll2007l) . 
An exact calculation of electrical conductivities, which are 
necessary for reconstructing the evolution of the magnetic 
field, would need to include the interaction between all 
the different particles. This manuscripts consider s only the 



electro n-ion contribution calculated by the code of lPotekhin 



(120081) . The values of the magnetic diffusivity and the Hall 
parameter used in this manuscript should be taken as order- 
of-magnitude estimates (see Figs. |9] and [TO]i. 
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Fig. 9 Microphysical magnetic diffusivity as a function of 
density for a temperature of 10^" K. 
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Fig. 10 The Hall parameter as a function of density for a 
temperature of lO^'' K. 

According to Fig. [TO] the Hall parameter and, therefore, 
the Hall effect vanish in deeper layers of a PNS's envelope 
even for a moderate magnetic field of 10^^ G. Only the outer 
layers with densities in the range of p ^ 10^*^-10^^ gcm^'^ 
are characterised by the Rb closer to unity. 

If the flow is laminar, Rm and S are orders of magni- 
tude larger than the examined parameter space in the sta- 
bility map of Fig. [T] due to the low diffusivity (see Fig.|9]l. 
If, however, the flow becomes turbulent, the turbulent dif- 



fusivity is much larger (iNaso et al. 2008), leadin g to much 
lower Rm (^10"*, as noted bv^ Kondic et alj201 l l and which 
are doable numerically), but still above the threshold for the 
onset of the Hall-MRI. An interesting feature of turbulent 
flows is that S is close to unity for densities larger than or 
equal to 10^^ g cm^^. These regions will be purely SHI un- 
stable, because the MRI is inhibited by diffusion. 

Therefore, the Hall-MRI instabilities shall set in after 
the seed field of the progenitor is sufficiently amplified. 
They would probably first destabilise the outer regions, per- 
haps spreading inward provided that the fields in the inte- 
rior become sufficiently strong (because the Hall parameter 
Rb scales linearly with the magnetic field). If the flow is 
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already turbulent for other reasons, even the pure SHI may 
be encountered in regions where the Lundquist number is 
small enough. 

Given that we have shown the growth rates of the in- 
stability to depend on a general orientation of the magnetic 
field, those multipolar field components which have the op- 
posite directions in different hemispheres, e.g. a quadrupole 
or a tilted dipole, will experience different growth. The 
fields with dominant negative orientation would be strongly 
restricted in their growth. This may have interesting impli- 
cations for observable properties of the field, with one hemi- 
sphere having different field strength and patterns than the 
other. 

The nonlinear effects related to this process are the self- 
coupling of the magnetic field through the quadratic Hall 
term, as well as the nonlinear nature of the flow/field inter- 
action where Maxwell stresses redistribute the angular mo- 
mentum in the fluid, braking the differential rotation on the 
time-scale comparable to the Alfven time (shorter than any 
time-scale related to viscosity effects). These effects will 
need to be taken into account in order to obtain the estimates 
of the final amplitude and topology of the magnetic field 
available through the Hall-MRI. The proximity of the ex- 
citation thresholds of axisymmetric and non-axisymmetric 
modes is another interesting outcome of this paper; the non- 
linear coupling of the modes may lead to different classes 
of surface field topologies with implications for the observ- 
ability of neutron stars. 
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